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We study a spin-1/2 system with Heisenberg plus ring exchanges on a four-leg triangular ladder using the den- 
sity matrix renormalization group and Gutzwiller variational wave functions. Near an isotropic lattice regime, 
for moderate to large ring exchanges we find a spin Bose-metal phase with a spinon Fermi sea consisting of three 
partially filled bands. Going away from the triangular towards the square lattice regime, we find a staggered 
dimer phase with dimers in the transverse direction, while for small ring exchanges the system is in a featureless 
rung phase. We also discuss parent states and a possible phase diagram in two dimensions. 



In a wide class of crystalline organic Mott insulators it is 
possible to tune from the strongly correlated insulating state 
into a metallic state. At ambient pressure such "weak Mott in- 
sulators" are perched in close proximity to the metal-insulator 
transition. The residual electronic spin degrees of freedom 
constitute a novel quantum system and can exhibit a myr- 
iad of behaviors such as antiferromagnetic (AF) ordering or 
a valence bond solid (VBS). Particularly exciting is the pos- 
sibility that the significant charge fluctuations in a weak Mott 
insulator frustrate the magnetic or other ordering tendencies, 
resulting in a quantum spin liquid. This appears to be re- 
alized in two organic materials [1-7] k-(ET)2Cu2(CN)3 and 
EtMe3Sb[Pd(dmit)2]2, both quasi-two-dimensional (2D) and 
consisting of stacked triangular lattices. Thermodynamic, 
transport, and spectroscopic experiments point towards the 
presence of many gapless excitations in the spin-liquid phase 
of these materials. 

The triangular lattice Hubbard model [1, 8-10] is com- 
monly used to describe these materials. At half filling the 
Mott metal-insulator transition can be tuned by varying the 
single dimensionless parameter, the ratio of the on-site Hub- 
bard U to the hopping strength t. On the insulator side at 
intermediate U/t, the Heisenberg spin model should be aug- 
mented by multispin interactions [11-16], such as four-site 
ring exchanges (see Fig. 1), which mimic the virtual charge 
fluctuations. Accessing a putative gapless spin liquid in 2D in 
such models poses a theoretical challenge. 

Slave particle approaches provide one construction of gap- 
less spin liquids and predict spin correlations that decay as 
power laws in space, oscillating at particular wave vectors. In 
the so-called "algebraic spin liquids" [17-20] these wave vec- 
tors are limited to a finite discrete set, often at high symmetry 
points in the Brillouin zone. However, the singularities can 
also occur along surfaces in momentum space, as they do in a 
"spinon Fermi sea" spin liquid speculated for the organic ma- 
terials [12-14]. We will call such a phase a "spin Bose-metal" 
(SBM) state [21, 22] to emphasize that it has metal-like prop- 
erties for spin and energy transport while the spin model is 
bosonic in character. 

It should be possible to access an SBM phase by system- 




FIG. 1: (color online). Picture of the Heisenberg plus ring Hamilto- 
nian on the four-leg ladder showing different two-spin and four-spin 
couplings. The isotropic model is defined by = Jg = J x + y = J, 
K x ,y = Kx,i+y ~ Ky^+y = K. We also study a broader 
phase diagram interpolating between the triangular and square limits 
by decreasing Jx+y [with appropriate scaling of the ring couplings 
Kx,x+ y = Ky,x+ y = (Jx+y / Jx) K x , y ]. The ladder has periodic 
boundary conditions in both directions. 

atically approaching 2D from a sequence of quasi- ID ladder 
models [21-23]. On a ladder the quantized transverse mo- 
menta cut through the 2D surface, leading to a quasi- ID de- 
scendant state with a set of low-energy modes whose number 
grows with the number of legs. These quasi-lD descendant 
states can be analyzed in a controlled fashion using numerical 
and analytical approaches. 

Heisenberg plus ring on a four-leg triangular ladder. — 
Pursuing this idea, we consider a spin-1/2 system with Heisen- 
berg and four-site ring exchanges, 

H = J2 2J ijSi-S j + tfp (7*1234 + H.c). (1) 

rhombi 

An earlier exact diagonalization (ED) work [24] on the 
isotropic 2D triangular lattice found that K > 0.1J destroys 
the 120° AF order. A subsequent variational study [12] sug- 
gested the spin Bose-metal phase for moderate to large K. 
A recent work pursued this model on a two-leg zigzag lad- 
der [22, 25] combining density matrix renormalization group 
(DMRG), variational Monte Carlo (VMC), and Bosonization 
approaches, and argued that it realizes a quasi- ID descendant 
of the SBM phase: a remarkable ID quantum phase with three 
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FIG. 2: (color online). Phase diagram of the Heisenberg plus ring 
model on the four-leg ladder interpolating between the triangular 
Jx = Jy = Jx+y = 1 and square J# = Jg = 1, Ji+ y = limits. 
The horizontal axis is the ring coupling K±. y while the vertical axis 
is the diagonal coupling J x + y , cf. Fig. 1; the other ring couplings are 
obtained according to Eq. (2). 



gapless modes and power law spin correlations at incommen- 
surate wave vectors that are the fingerprints of the parent 2D 
phase. 

The two-leg ladder is still far from 2D. We take a signifi- 
cant step and study the model on a four-leg ladder. We first 
consider the case where all nearest neighbor bonds have the 
same coupling J and all rhombi have the same coupling K; 
thus there is a single parameter K/J. 

We study the model numerically using DMRG/ED com- 
bined with VMC calculations. All calculations use periodic 
boundary conditions. The DMRG calculations keep m = 
3600-5000 states per block [26-28] to ensure accurate results, 
and the density matrix truncation error for our systems is of 
the order of 10 -5 (typical relative error for the ground-state 
energy is 10~ 3 or smaller). Information about the state is ob- 
tained by measuring spin, dimer, and (scalar) chirality struc- 
ture factors. 

The phase diagram from such a study using 12x4 and 18x4 
ladders can be seen in Fig. 2; the isotropic case is the horizon- 
tal cut at Jx+y 1 1 Jx = 1. For small K/J < 0.15 the system 
is in a rung phase, whose caricature can be obtained by al- 
lowing Jy ^S> J x , J x +y where the rungs effectively decouple. 
This phase is gapped and has only short-range correlations. 
In the model with isotropic couplings the rungs have rather 
strong connections: we find that the x and x + y bonds have 
more negative Heisenberg energies than the y bonds. Nev- 
ertheless, the data suggest that the system is in a featureless 
gapped phase. A further test is provided by increasing J y from 
the isotropic case, and we indeed observe a smooth evolution 
in all measurements towards the strong rung phase. 

Near K/J = 0.2-0.25, the DMRG ground state breaks 
translational symmetry. The pattern obtained on both the 
12 x 4 and 18 x 4 systems is illustrated in Fig. 3(a). This 
state has strong y bonds forming columns along the ladder 
direction and strong x + y bonds in the connecting arrange- 
ment. Note that we also expect a degenerate state depicted in 




FIG. 3: (a) Symmetry breaking pattern found in DMRG on the 
isotropic system at K/J — 0.2-0.25. (b) State degenerate with (a) 
in the presence of x -f-> x + y symmetry. Both (a) and (b) can be 
viewed as triangular VBS states with dimers on the y bonds but dif- 
ferent column orientations. The staggered patterns on the x + y and 
x bonds correspondingly are expected on the triangular lattice and 
follow a rule that each triangle contains only one strong bond. Upon 
going to the square limit by decreasing Jx+y, we find state (b), which 
connects to a staggered j/-dimer state. 

Fig. 3(b), since the x and x + y directions are equivalent on 
the isotropic ladder. The states shown in Fig. 3 are a subset of 
possible VBS states on the isotropic 2D triangular lattice, and 
the selection must be due to the finite transverse size. The se- 
lection of (a) in the DMRG must be due to symmetry breaking 
terms that exist in the way it is building up the multileg sys- 
tem. Such terms are tiny and translationally invariant ground 
states are obtained for all other phases without intrinsic de- 
generacy (we also verified that the DMRG obtained identical 
results to the ED for 8 x 4 systems). 

SBM phase. — For K/J > 0.3, we do not find any pattern of 
bond ordering in real space and no indication of Bragg peaks 
in the dimer or chiral structure factors. The correlation func- 
tions are also markedly different from the rung phase at small 
K. The 12 x 4 and 18 x 4 systems remain in essentially the 
same state for a range of control parameters 0.3 < K/J < 1. 
Thus, a putative spin-liquid phase is established based on 
finite-size analysis of the DMRG results. Spin and dimer cor- 
relations are rather extended in real space and show complex 
oscillations. The momentum space structure factors allow a 
more organized view and show many features that, remark- 
ably, can be manifested by simple variational wave functions 
for the SBM phase. 

To this end, we perform a VMC study using spin-singlet 
trial wave functions that can be viewed as Gutzwiller projec- 
tions of spinon hopping mean field states. More systemati- 
cally, we vary directly the shape of the "spinon Fermi sea" 
in the momentum k space. There are four transverse val- 
ues k y — 0, ±7r/2, 7t, and for each we can allow an arbitrary 
"Fermi segment," i.e., a contiguous region of occupied fc or- 
bitals. For the 8x4 system, we optimized the trial energy over 
all distinct locations of these segments, the only restriction be- 
ing the specified total filling, and found that only three bands 
are populated in a manner that respects the lattice symmetries. 
For the 18 x 4 system, from the outset we restricted the op- 
timization to such three-band states with inversion symmetry 
and found a state shown in Fig. 4. 

Figure 5(a) shows the spin structure factor measured in 
DMRG and calculated using the optimal VMC state on the 
18 x 4 ladder, while Fig. 5(b) shows the dimer structure fac- 
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FIG. 4: "Spinon Fermi sea" of the optimal VMC state for _ftT/ J > 
0.3 on the 18 x 4 ladder. Circles denote single-particle orbitals in 
fe space; f and J, spinons occupy orbitals shown with filled circles. 
There are three partially filled bands and we refer to the resulting 
state as SBM-3. The dotted line indicates k x = 0, where we are 
implicitly using antiperiodic boundary conditions in the horizontal 
direction. This choice, as well as the location of k x — 0, is actually 
arbitrary as the Gutzwiller-projected wave function is invariant under 
a global shift of the Fermi sea. 



tor for bonds oriented in the y direction. In both figures, 
we see sharp peaks at wave vectors (3 x 27r/18,7r/2) and 
(10 x 2tt/18, 7t/2). These are reproduced by the VMC wave 
function and the wave vectors can be associated with spinon 
transfers between right-mover and left-mover Fermi points 
differing by Ak y — ±tt/2 in Fig. 4. 

Figure 5(c) shows the dimer structure factor for bonds ori- 
ented in the x + y direction. Here, we see notable peaks at 
wave vectors ±(2 x 2-7r/18,7r) and ±(6 x 27r/18,7r) in ad- 
dition to some of those previously mentioned. The VMC 
agreement at these two points is again striking and can be as- 
sociated with ±2kp 7r/ ' 2 ^ and ±2kp~ n ^ spinon transfers in 
Fig. 4. While there are quantitative discrepancies between the 
VMC and DMRG approaches, the overall agreement in the 
location of sharp features is notable. We do not show chiral- 
ity structure factor results but can report that no ordering was 
observed over the regime of parameters shown in our phase 
diagram. 

A summary of the gauge field theory for the SBM-3 phase 
is as follows. In the mean field, there are three partially filled 
bands for each spinon species as in Fig. 4. Linearizing near 
the Fermi points and bosonizing, there are total of six modes: 
three in the "charge" sector and three in the "spin" sector. In 
the quasi- ID system, the gauge fluctuations beyond the mean 
field eliminate the overall charge mode, leaving five gapless 
modes. In principle, one should be able to verify this predic- 
tion by analyzing the entanglement entropy with DMRG and 
extracting the central charge. Indeed, we were able to confirm 
this for an 8 x 4 system, but exhausted the DMRG resources 
for larger sizes. 

The three modes in the spin sector cannot have nontriv- 
ial Luttinger parameters due to the SU(2) invariance, while 
there are two Luttinger parameters describing the charge sec- 
tor. Inspection of all allowed interactions shows that the SBM- 
3 can in principle be a stable phase, although there are many 
channels where this multimode system can become unstable. 
Dominant spin and dimer correlations are expected at wave 
vectors ±2kp a , where kp a runs over the three right-mover 
Fermi points in Fig. 4. Potentially enhanced correlations are 
also expected at wave vectors ±(kp a + kpb)\ at each such 
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FIG. 5: (color online). DMRG and VMC structure factors for the 
18 x 4 triangular ladder at K/J = 0.6. The horizontal axis is a 
linear representation of the L x 4 mesh of (q x ,q y ) points, and the 
disjoint curves from left to right correspond to q y = 0, 7r/2, n, and 
3tt/2. (a) Spin structure factor, (b), (c) Dimer structure factor for 
bonds oriented in the y direction and x + y direction. 



wave vector, the "multiple!" of observables with the same 
power law contains spin, dimer, scalar chirality, and vector 
chirality. The Gutzwiller wave functions correspond to a fine- 
tuned theory where all Luttinger parameters are equal to unity 
and the above wave vectors show the same power law x~ 5 ^ 3 . 

While the DMRG does not show some of the expected wave 
vectors, the overall match with the VMC suggests that this 
may be due to matrix element effects. It can also be that the 
SBM-3 is eventually unstable here, but it clearly is a good 
starting point for understanding the remarkable phase found 
in the Heisenberg plus ring model. 

Interpolation between the triangular and square limits. — 
Motivated by relatives of the k-(ET) 2 Cu 2 (CN)3 and 
EtMe3Sb[Pd(dmit)2]2 spin-liquid materials, we extend the 
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study by allowing a different diagonal coupling J x + y < J x — 
J y . We also vary K XtV , while the remaining ring couplings 
are fixed by 

K x ,x+y Ky,x+y j -Kx,y • (2) 

In an anisotropic electronic system with hoppings t„ in the a 
direction, the Heisenberg couplings are Ja <~ t\/U and the 
ring couplings K & g ~ tltf/U 3 , so their "anisotropies" are 
indeed related: g ~ J a Jg/£7. 

Figure 2 gives the phase diagram in the K XtV /J x — 
Jx+y/Jx plane determined from the DMRG and VMC. We 
see three prominent phases. Along the K x ^ y = axis, the sys- 
tem is in the rung phase. Going to the square limit, J x + y = 0, 
and then increasing K XiV , the system undergoes a transition to 
a staggered dimer phase for K x ^ y > 0.8. This agrees with an 
earlier study [29] of the 2D square lattice model with ring ex- 
changes that found staggered dimer phase in the same regime. 
In the present four-leg system, the dimers orient transverse to 
the ladder (y direction). As J x + y increases towards the tri- 
angular regime, the dimer phase expands to smaller values of 

Kx,y' 

Significantly, the dimer phase disappears for anisotropy 
0.8 < J x+y < 1 and moderate to large ring exchange. Here 
the DMRG finds a spin-liquid state that fits with the SBM-3. 
The dimer phase touches the triangular axis near K = 0.2- 
0.25 where it becomes degenerate with the VBS state found 
earlier, cf. Figs. 3(a) and 3(b). 

Discussion. — The four-leg ladder captures a good deal of 
local physics of the 2D model and allows guesses about the 2D 
phase diagram. First, in the rung phase near the square limit 
we observe strong spin correlations at (it, it) and can view this 
region as a ladder descendant of the square lattice Neel state. 
Second, the staggered y dimer state is a descendant of the 2D 
staggered dimer state. In 2D there is a degeneracy between 
cases where the strongest bonds are oriented in the y or x di- 
rection. Third, the four-leg SBM-3 state is a direct descendant 
of the 2D SBM. To summarize, we expect the 2D Neel, stag- 
gered VBS, and spin-liquid phases to occupy roughly similar 
regions as the rung, staggered y dimer, and SBM-3 phases in 
Fig. 2. We do not venture to speculate how the three phases 
meet, particularly since additional phases enter into competi- 
tion. Specifically, the 2D triangular lattice with K < 0.1 has 
the 120° AF phase [24]. Also, series expansions [30] suggest 
that the Heisenberg model (K = axis) has a columnar VBS 
phase for 0.7 < J x + y /J x < 0.9, which is different from the 
staggered VBS stabilized by the ring exchanges. 

Naive Hubbard model estimates for the organic spin-liquid 
materials suggest that they lie in the challenging regime of 
strong frustration (J x + y ~ J x ~ J y ) and small to interme- 
diate K ~ 0.2J. Our study shows that the SBM phase is 



a viable contender. More realistic treatments such as inclu- 
sion of further ring exchanges (pursued systematically for the 
Hubbard model in [16]) and long-range Coulomb interactions 
may tilt the balance towards the spin-liquid phase and deserve 
further study. 
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